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Abstract 

We consider two variations of the discrete car parking problem where at every 
vertex of Z a car arrives with rate one, now allowing for parking in two lines. 

a) The car parks in the first line whenever the vertex and all of its nearest neighbors 
are not occupied yet. It can reach the first line if it is not obstructed by cars already 
parked in the second line ( "screening" ) . 

b) The car parks according to the same rules, but parking in the first line can not 
be obstructed by parked cars in the second line ("no screening"). In both models, a 
car that can not park in the first line will attempt to park in the second line. If it is 
obstructed in the second line as well, the attempt is discarded. 

We show that both models are solvable in terms of finite-dimensional ODEs. We 
compare numerically the limits of first and second line densities, with time going to 
infinity. While it is not surprising that model a) exhibits an increase of the density 
in the second line from the first line, more remarkably this is also true for model b), 
albeit in a less pronounced way. 

AMS 2000 subject classification: 82C22, 82C23. 

Key-Words: Car parking problem, Random sequential adsorption, Sequential fre- 
quency assignment process, Particle systems. 

1 Introduction 

Car parking, first considered in a mathematical way by Renyi jlOj in 1958, gives rise to 
interesting models that in several variations have been applied in many fields of science. 
In the original car parking problem, unit length cars are appearing with constant rate in 
time and with constant density in space on the line where they try to park. A new car is 
allowed to park only in case there is no intersection with previously parked cars. Otherwise 
the attempt is rejected. Renyi proved that the density of cars has the limit 0.7475..., the 
so-called parking constant. In the simplest discrete version of the car parking problem, 
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cars of length 2 try to park at their midpoints randomly on Z. This model has been solved 
analytically as well [3]. 

This model belongs to a wider class of more complicated models of deposition with 
exclusion interaction. Usually such models are not analytically solvable. In physical 
chemistry "cars" become particles which are deposited in layers on a substrate, a process 
called random sequential adsorption (RSA). A variety of related models are studied. For 
a review of recent developments see . Moreover, models with more complicated graphs 
e.g. (random) trees have been investigated [2 ,[9j,[6],[5]. 

Multilayer variations of the model are used to describe multilayer adsorption of parti- 
cles on a substrate [8] and the sequential frequency assignment process [3 which appears 
in telecommunication. In these papers it is also observed that the density in higher layers 
increases up from the first layer, which at first seems rather counterintuitive. Heuristic 
arguments for monotonicity of densities were found in specific models [8!, but no rigor- 
ous proofs could be given yet. Moreover Privman finds numerically a scaling behavior 
of the density in a similar RSA model [8^ with slightly different adhesion rules which is 
notoriously difficult to explain mathematically. 

In the present paper we aim for a rigorous investigation and treat two versions of the 
discrete two- line car-parking problem with cars of length 2. First we describe the dynamics 
of the car parking process without screening and also with screening. Then we provide 
the solutions of these models by reducing them to closed finite dimensional systems of 
ODEs for densities of local patterns, see Theorems [T] and [2j That it is possible to find 
a finite-dimensional dynamical description is quite remarkable. It is not obvious, and in 
fact our method ceases to work for a three-line extension of the model without screening 
where an infinite system appears. 

A second remarkable fact is that, even without screening, the second line density is 
higher than the first. Cars do not communicate or plan a common strategy and their 
arrival is random, but they seem to use the resources in the second line more efficiently, 
once they have been rejected in the first line. 



2 The Dynamics 

We will define a Markov jump process on the (suitably coded) occupation numbers m = 
(mi)iez gO = {0,1,2,3}^. 

Here the spin rui denotes the joint occupation numbers at vertex i at height 1 and 2. 
It is useful for short notation to interpret the occupation numbers at various heights as 
binary digits and write ordinary natural numbers. That is we write 

if vertex i is vacant in the first and second line 

1 if vertex i is occupied in the first but not in the second line 
* 2 if vertex i is occupied in the second but not in the first line 

3 if vertex i is occupied in the first and in the second line 

so that rrii E {0, 1,2, 3}. The dynamics of the process is defined in terms of the generator 
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which is given by the right hand side of the differential equation 
Jt 



E™/(m(t)) = [fim'''') - f{m) c(s; mfc_i, m^, mfc+i) 



(2.1) 



with 



s,k 



rrii \i k ^ i 
s ii k = i 



denoting the configuration which has been obtained by m by changing the configuration 
in i to s. Here E™ denotes the expected value with respect to the process, started at the 
initial configuration m. 

Two-line parking rates 

The rates are either equal to zero or one. They are 1 precisely in the following cases. 

1. 1-^ 1 Adding a car in the first line at site i. For the model without screening we 
have 



c(l;0,0,0) = c(l;2,0,0) = c(l;0,0,2) = c(l;2,0,2) = 1 



(2.2) 



Indeed, this occurs when the site itself is empty on the first and second line and the 
nearest neighbors are empty in the first line, see figure [T] for an example. 



i 



-1 1 



-1 1 



-1 1 



-1 1 



Figure 1: Configurations of vertices —1,0, and 1 that allow a transition from mo = to 
mo = 1 in the model without screening. In the model with screening only the most left 
configuration allows a transition to niQ = 1. 

In the screening model however, cars in the second line will obstruct cars from reach- 
ing the first line. Therefore, in the screening model we have as the only nonvanishing 
rate 

cS^(l;0,0,0) = 1 (2.3) 

2. 1-^ 2 Adding a car in the second line at i while the first line was empty at the site 

c(2; 1,0, 0) = c(2; 0, 0, 1) = c(2; 1, 0, 1) = 1 (2.4) 

Indeed, this occurs when there was a supporting site z — 1 or i + 1 or both with one 
car in the first line. This is true for both models. 
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3. 1 3 Adding a car in the second line while the first line was full at the site 



c(2;0,l,0) = l (2.5) 

Indeed, this occurs when there are no obstructing cars right and left at height 2. 
There can be no obstructing cars right and left at height 1 because there could not 
be a car in the first line at i otherwise. This is true for both models. 

All other transitions are impossible. 

This generator defines a Markov jump process on the infinite graph by standard theory 



[T], such that (2.1 1 holds for any local function / : Q 



3 Results 

We provide a closed system of differential equations for the densities of occupied sites, 
involving densities of finitely many local patterns, in both models. First we need some 
definitions. Here and in the following we use for the densities at single sites, and triples 
of sites the notation 

Dt{s) := Pt{mo = s) 
Dt{s-i,so,si) := Pt{m-i = S-i,mo = so,mi = si) (3.1) 



Further we need the following "one-sided densities" 

/t(s) :=Pt(mi = s|iVo(i) = 0), for s = 0,1,2 
Rt := Pt{mi = 1, ma = 0|iVo(0 = 0) 

where Nj{t) denotes the Poisson counting process of events of car arrivals at site j. 

As our main result we show that the time-evolution of these densities gives rise to a 
closed ODE. 



Theorem 1 Two-line Parking without Screening. 

The time evolution of the probability vector {Dt{0), Dti^), Dt{2), Dt{3)) obeys the following 
system of differential equations. 

|a(0) = -(/*(0) + ft{2))'e-' - (2/,(0)/i(l) + ft{lf)e-' 
^ A(l) = iftiO) + /t(2))2e-* - A(0, 1,0) 

(3.3) 

|A(2) = (2/t(0)A(l) + /i(l)2)e-* 
^A(3) = A(0,1,0) 
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with initial conditions Dq{s) = ls=o, where the vector {ft{0),ft{l),ft{'^),Rt) obeys the 
linear ODE 



d_ 

dt 
d_ 

'dt 
d_ 

dt 



MO) = -Mo)e~' - Miy 



m = me~' + me 



-t 



Rt 



m = ft{l)e 



-t 



d_ 
dt 



Rt = /t(0)(e-* - - f,{l)te- 



■2t 



Rt 



with initial conditions (/o(0), /o(l), /o(2), i?o) = (1,0,0,0), 
and finally, -Di(0,l,0) is obeying the equation 



d_ 
dt 



A(0, 1, 0) = ftiOfe-' - A(0, 1,0) - 2RtM0)e-' - 2Rtft{l)e-' 



with Do(0,l,0) = 0. 



(3.4) 



(3.5) 



The system above can be solved numerically and the results are depicted in figure [2} 
As it can be seen in the right figure, surprisingly the value of Dt{2) has a slightly higher 
limit than -Dt(l). This clearly means that the second line has a higher limit density of cars 
than the first line. This result is independently confirmed by simulations of the parking 
process measuring the empirical densities. 



Densities on the First and Second Line 
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.a- 

i 0.6 



- - D,(0) 
■ D,(1) 

D,(2) 

D,(3) 



0.3^ 
0.295 - 

0.29 - 
0.285 - 

0.28 - 
0.275 - 

0.27 - 
0.265 - 

0.26 - 

0.255 - 

0.25 — 




D,(0) 
D,(1) 
- D,(2) 



Figure 2: Numerical solution of the densities in the first and second line. The right figure 
zooms in on the limits of Dt(l) and Dt(2) to show they tend to different values. 

A similar system of equations holds for the model with screening. Recall that in this 
model cars are not allowed to pass cars on the second line to reach a void on the first line. 
This results in less possibilities of filling voids of the first line than in the model treated 
above. In fact we can derive the ODEs of this model by simply deleting those terms in 
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(3.4) that represent the possibihty of "jumping" over a car in the second hne to reach a 



void on the first hne. So, we get 

Theorem 2 Two-line Parking with Screening 

The time evolution of the probability vector (I?|'^(0), Z)|'^(2), L>|'^(3)) obeys the fol- 

lowing system of differential equations 

|z^r(o) = -(/r(o) + /r(i))v* 

^A^^(i) = /r(o)V*-Ami,o) 

^ (3.6) 

-Df{2) = (2/r(o)/r(i) + /r(i)')e-* 
j^Dr{3) = Dr {0,1,0) 

with initial conditions Dq^(s) = ls=0; where the vector {f^^(0), f^^{l), f^'^{2), Rf^) obeys 
the linear ODE 

j^fm = -fme-'-fme-' 



|/r(i) = /r(o)e-*-i?f 
|/r(2) = /r(i)e-* 



(3.7) 



|i?r = /r(0)(e-* - te-2*) - /r(l)te-2* - Rt 



with initial conditions (fg^O), f§%l), f§%2), R'q'') = (1,0,0,0), 
and finally, -D|^(0,1,0) is obeying the equation 

jDnO, 1,0) = fr{Ofe-' - DfiO, 1,0) - 2i?f /r(0)e-* - 2i?f /r(l)e-* (3.8) 
with D'q%0,1,0) = 0. 

4 Proofs of Theorem 1 and Theorem 2 

The following lemmas are used to prove our theorems. 

Lemma 4.1 The probability vector {Dt{0), Dt{l), Dt{2), Dt{3)) obeys 

^ A(0) = - A(0, 0, 0) - 2A(2, 0, 0) - A(2, 0, 2) - 2A(1, 0, 0) - A(l, 0, 1) (4.1) 

^A(l) = A(0,0,0)+2A(2,0,0)+A(2,0,2)- A(0,1,0) (4.2) 

^A(2) = 2A(1,0,0)+A(1,0,1) (4.3) 

|a(3) = A(0,1,0) (4.4) 
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Remark: Summing over the four right hand sides we get zero, due to the fact that 
we have summed a probabihty vector. It is also interesting to check that 

^ A(l) + ^ A(3) = A(0, 0, 0) + 2A(2, 0, 0) + A(2, 0, 2) (4.5) 
recovers the ODE for the density in the first hne. 

Proof: Fix an arbitrary vertex. Let us call this vertex 0. Starting from the dynamics 
(2.1) and using symmetries we have 

^ A(0) = -A(0, 0, 0) - 2A(2, 0, 0) - A(2, 0, 2) - 2A(1, 0, 0) - A(l, 0, 1) (4.6) 

Indeed, the first three terms correspond to adding a car in the first line, the next two 
terms correspond to adding a car in the second line, see figure |3j 



A(0,0,0) A(2,0,0) A(0,0,2) A(2,0,2) A(1,0,0) A(0,0,1) A(1,0,1) 



Figure 3: List of all occupancy configurations of vertices -1, and 1 in the first and second 
line that (may) contribute to a reduction of the proportion of Dt{0). 

The other three differential equations are derived in a similar way. 
Lemma 4.2 The triple site densities Dt{s,0, s') and the one-sided densities as defined in 



3. 1 and \3.^ respectively are related in the following way 

Dtis,0,s') = ftis)ftis')e-' (4.7) 

for{s,s') G {(0,0), (0,1), (1,0), (1,1), (2,2)}. 

Proof: We note that for the mentioned choices of (s, s') conditioning on non-arrival 
at zero does not change the probability, that is 



Dtis, 0, s') = Pt{m^i =s,mo = 0, mi = s' , No{t) = 0) 
= i-f(m„i = s,mo = 0,mi = s'\NQ{t) = 0)e~ 



(4.8) 



In the next step we note that, conditional on the event that no car has arrived at the site 
0, the dynamics for the two sides that are emerging from is independent. Consequently 
we have 

Pt{m-i = s,mo = 0,mi = s'\No{t) = 0) 

= Pt{mi = s\No{t) = 0)Pt{mi = s'\No{t) = 0) (4.9) 
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This concludes the proof of the Lemma. □ 
Next we look at the time-evolution of the "one-sided densities" . 

Lemma 4.3 The vector {ft{0),ft{l),ft{2),Rt) obeys the ODE 
jjm = -/t(0)e-* - /t(l)e"* - /i(2)e-* 



= me-' + ft{2)e-' - Rt 



d 
dt 

|/,(2) = me-' 



(4.10) 



jRt = hme-' - - mte-^' " Rt 

with initial conditions ft^Q^s) = ls=o and Rt=o = 0. 

Remark 1: Note that combining the equations of ^ft{0) and ^ft{2) readily gives 

A(0) + /t(2) = exp(e-* - 1) (4.11) 

which is a known result for the first line in a semi-infinite chain [31. 

Remark 2: Note also that because Rt := Pt{mi = 1,^2 = 0|A'o(t) = 0) we have in 
fact Rt = f^m. 

Proof: To derive ODEs for these densities we employ the generator of the process, 
while putting the term at the site to sleep, and correspondingly the spin at zero to be 
the constant ttt-q = 0. For the first quantity we get 

|/,(0) = |p,(mi = 0|iVo(t) = 0) 

= -Pt{mi = 0,m2 = 0\No{t) = 0) - Pt{mi = 0,m2 = 2\No{t) = 0) 
-Ptimi = 0,m2 = l|iVo(t) = 0) 
= -Pt(mi = 0,m2 = 0\No{t) = 0,Ni{t) = 0)6"* 

- Pt{mi = 0,m2 = 2\No{t) = 0, Ni{t) = 0)e-' (4.12) 
-Pt{mi = l|iVo(t) = 0)e-* 

= -Ptimi = 0\No{t) = 0)e-* - Pt{mi = 2\No{t) = 0)6"* 

-Pt{mi = l|iVo(t) = 0)e"* 

= -ft{0)e-' - ft{l)e-' - me-' 



Next we have 



|/,(l) = |p,(mi = l|7Vo(t) = 0) 

= Pt{mi = 0, m2 = 0|7Vo(i) = 0) + Pt{mi = 0, m2 = 2|iVo(t) = 0) 

-Pt(mi = l,m2 = 0|7Vo(t) =0) ^^-^^^ 

= /t(0)e-* + /t(2)e-* - Rt 



8 



and 



|/i(2) = |p,(mi = 2\No{t) = 0) 



Ptimi = 0, ms = l|iVo(i) = 0) (^"^^^ 
= Pt(mi = l|iVo(i) = 0)e-* 

Finally, we get 

= ^^Pt{mi = l,m2 = 0\No{t) = 0) 

= Pt{mi = 0,m2 = 0\No{t) = 0) - Pt{mi = l,m2 = 0\No{t) = 0) (4.15) 

- Pt{mi = l,m2 = 0, ma = 0|A^o(i) = 0) 

- Pt{mi = l,m2 = 0,m3 = l\No{t) = 0) 

Using conditioning on non-arrival again we get 
^Pt(mi = l,m2 = 0\Noit) = 0) 

= Pt{mi = 0\No{t) = 0)e-* - Pt{mi = 1, ms = 0\No{t) = 0) (416) 

- Pt{mi = l\No{t) = 0,N2{t) = 0)Pt{mi = 0\No{t) = 0)6"* 

- Ptimi = l\Noit) = 0,N2{t) = 0)Pt(mi = l\No{t) = 0)6"* 

Clearly we have 

Ptimi = l\No{t) = 0,N2{t) = 0) = te-' (4.17) 

because there is precisely one car at 1 if and only if precisely one car arrived conditioning 
on no cars at and 2. This shows that the last ODE is correct and concludes the proof 
of the lemma. □ 
The only remaining term whose time-evolution we need to consider is Dt(0, 1,0). 

Lemma 4.4 D((0, 1,0) is a solution of the differential equation 

|a(0, 1,0) = ftiOfe-' - A(0, 1,0) - 2Rtft{0)e-' - 2Rtft{l)e-' (4.18) 

Proof: We note that 

^ A(0, 1,0) = A(0, 0, 0) - A(0, 1,0) - 2A(0, 1, 0, 0) - 2A(0, 1,0, 1) 

The first term is for adding a car at the central site from the vacuum, the second for 
adding a car at the central site at height one. The last two terms are for adding a car to 
the right of the central site. As we already know we have 

A(0, 0, 0) = Pt{mi = 0\NQ{t) = Ofe-^ (4.20) 
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Using conditioning on non- arrival at 2 we get, by reflection invariance 

A(0, 1, 0, 0) = Pt(mi = 1, ma = 0| A^oW = 0)Pt{mi = 0\No{t) = 0)6"* 
= RtMO)e-' 

For the last term we get in the same way 

A(0, 1, 0, 1) = Pt{mi = 1, ma = 0\No{t) = 0)Pt{mi = l\No{t) = 0)6"* 



(4.21) 



(4.22) 



□ 



Proof of Theorem [TJ Combining the results of lemma 4.1 4.2, 4.3 and 4.4 proves 
Theorem 1. □ 



Proof of Theorem [2j The proof follows from Theorem [T] by deleting every term 
that represents the possibility of skipping a second line car to reach a void in the flrst line. 



This results in deleting /t(2)e from the flrst two equations of 4.3 and Dt{2,0,0) and 



A(2,0,2) from 4.1 andO 



□ 



Density of the first line 



Density of the second line 




- density without screening 
density with screening 




- density without screening 

- density with screening 
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Figure 4: Densities of the first and second line. 



5 Conclusion 

We introduced two extensions of the classic parking problem to a two-line model i.e. a 
model with screening and a model without screening. For both models we derived closed 
systems of finite-dimensional ODEs from which the time-evolution of the densities in the 
first and second line can be obtained. Interestingly, the numerical solution of the ODE 
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shows that the final densities in the second line are higher than those in the first line, for 
both models. The increase factor in the model without screening is approximately 



t^oo Dt{l) + A(3) 0.432332 

It is known by analytical computations ^ that -Dt(l) + Dt{3) approaches (1 — e^^)/2 ~ 
0.432332 for large t, which provides a checkup for the numerics. In the screening model 
we find 

i^oo Lff (1) + L>f (3) 0.366475 

In other words, in both models the cars seem to exploit the resources in the second 
line in a (slightly) more efficient way than in the first line. 
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